Atomic scale elastic textures coupled to electrons in superconductors 
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We present an atomic scale theory of lattice distortions using strain related variables and their 
constraint equations. Our approach connects constrained atomic length scale variations to continuum 
elasticity and describes elasticity at all length scales. We apply the general approach to a two- 
dimensional square lattice with a monatomic basis, and find the atomic scale elastic textures around 
a structural domain wall and a single defect, as exemplary textures. We clarify the microscopic origin 
of gradient terms, some of which are included phenomenologically in Landau- Ginzburg theory. The 
obtained elastic textures are used to investigate the effects of elasticity-driven lattice deformation on 
the nanoscale electronic structure in superconductor by solving the Bogliubov-de Gennes equations 
with the electronic degrees of freedom coupled to the lattice ones. It is shown that the order 
parameter is depressed in the regions where the lattice deformation takes place. The calculated 
local density of states suggests the electronic structure is strongly modulated as a response to 
the lattice deformation — the elasticity propagates the electronic response over long distances. In 
particular, it is possible for the trapping of low-lying quasiparticle states around the defects. These 
predictions could be directly tested by STM experiments in superconducting materials. 

PACS numbers: 81.30.-t, 74.25. Jb, 74.50.+r, 61.50.Ah 



In many complex electronic materials such as cuprates, 
manganites, ferroelastic martensites, and titanates, un- 
expected multiscale modulations of charge, spin, polar- 
ization, and strain variables have been revealed by high 
resolution microscopy 1]. It is increasingly evident that 
the nonuniform textures found in these doped materi- 
als have intrinsic origins: they arise from coupling be- 
tween various degrees of freedom. The textures funda- 
mentally affect local and mesoscopic electronic, magnetic 
and structural properties, which are central to the func- 
tionality of correlated electronic materials. There is am- 
ple evidence for significant coupling amongst the elec- 
tronic degrees of freedom with the lattice distortions in 
cuprates, manganites, and ferroelectrics. The charge car- 
rier doping can act as a local stress to deform surround- 
ing unit cells 1]. We might employ a Landau-Ginzburg 
(LG) theory to study the coupling between the electronic 
(Cooper pair) and lattice (strain tensor) degrees of free- 
dom in superconductors. However, the LG theory can 
only rigorously describe the long wavelength behavior. 
New generations of experimental tools to probe individ- 
ual atoms and local environments 2] and the aforemen- 
tioned growing interest in complex functional materials, 
emphasize the importance of accurately describing the 
local electronic properties and lattice distortion at the 
atomic scale. In this work, we first present a microscopic 
description of elasticity. We introduce appropriate inter- 
cell and intra-cell distortion modes and show how the 
form of the elastic energy recovers the correct phonon 
spectra. The discreteness of the lattice, choice of modes 
and constraints among them give rise to an anisotropic 
gradient expansion for the elastic energy. This leads to 
interesting elastic domain wall and defect textures. We 
then couple these textures to the electronic degrees of 
freedom and study microscopically the influence of strain 
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FIG. 1: Normal distortion modes for a square object of four 
atoms in 2D. 



on electronic wavefunctions in both s and d-wave super- 
conductors. 

Our approach is general, but we illustrate it here in 
detail for the simplest case, namely a square lattice in 
two-dimensional (2D) space with a monatomic basis. We 
find that the most convenient strain-related variables for 
atomic scale distortions are the normal distortion modes 
of an elementary square object of four atoms (Fig. 
The first three distortion modes in Fig. Q correspond to 
the usual dilatation (ei), shear (62), and deviatoric (es) 
strains of the continuum elasticity theory for a square 
lattice 3j . The next two degenerate modes in Fig. ^ 5+ 
and 5_, correspond to the "intracell" or "shuffle" modes 
of the square lattice 1^, which are absent in continuum 
elasticity theory. Our approach uses these five distortion 
variables defined for each plaquette of four atoms at z, 
i + (10), i + (11), and i + (01), where i represents the 
coordinate of the lattice points, to describe the elastic 
energy 5]. 

Since the five variables are derived from two displace- 
ment variables for each lattice site, they are related by 
three constraint equations. By representing ei, 62, 63, 
s+, and S- in terms of displacement variables and 
in k (wavevector) space and eliminating and d^ ^ the 
constraint equations are obtained. One of them is the 
microscopic elastic compatibility equation, which relates 



2 



strain modes: 

(1 — cos kx cos ky)ei{k) — sin kx sin kye2{k) 
-\-{coskx — cosky)es{k) = 0. (1) 

The other two relate the intracell and the strain modes: 

2cos^cos^s±(^) T^sin (^t^^^ ^lik) 

±zsin(^^)e3(^) =0. (2) 

These constraints generate anisotropic interactions (from 
the lattice symmetry) between atomic scale strain fields, 
similar to the compatibility equations in continuum the- 
ory 3], but now including the intracell modes. In the 
long wavelength limit, our description naturally repro- 
duces the continuum results: For ^ ^ 0, the above con- 
straint equations can be written in real space as 

V'ei(r) - 2V,V^e2(r) + (Vg - Vl)e^{r) = 0, (3) 
s±{r) = \ [(V, ± Vx) ei(r) + (V, T V,) es{r)] .(4) 

Equation Q is the usual compatibility equation in con- 
tinuum theory. Equation shows that the spatial varia- 
tions of strains always generate intracell modes^ the mag- 
nitudes of which vanish as the inverse of the length scale 
of the strain mode variations. In continuum LG theory, 
the energy associated with the gradient of strains is re- 
sponsible for domain wall energies as, e.g., in structural 
phase transitions 4] . The above result shows that the in- 
tracell modes are the origin of such energy terms. Since 
our strain-related variables become identical to conven- 
tional strain variables for A: ^ , various length scale 
lattice distortions may be described in a single theoret- 
ical framework. This makes it possible to study typ- 
ical multiscale situations where both short- and long- 
wavelength distortions are important. It also provides 
a natural framework for incorporating interactions be- 
tween atomic scale strain-related fields coupled to other 
degrees of freedom in functional materials (below). 

The following analysis of the simple harmonic elastic 
energy for the square lattice further exemplifies the util- 
ity of these variables. We consider the simplest energy 
expression by approximating the total elastic energy as 
the sum of the elastic energy of each square: 

l n=l,2,3 m=+,- 

(5) 

where and B denote elastic moduli and 'intracell mod- 
ulus', respectively. Since some of the atomic pairs are 
shared by two square plaquettes of atoms, the parame- 
ters in Eq. (5) should be appropriately renormalized. A 
robust way to determine the parameters is to compare the 
phonon spectrum of our model with experimental data. 
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FIG. 2: An example of calculated phonon spectra (a) with 
and (b) without intracell modes for a 2D square lattice with 
a mono- atomic basis. The upper phonon branch is shown for 
both cases. (M — Ti — V). 

The phonon spectrum has been obtained • A typical 
spectrum (upper branch) for Ai=5, ^2=4, ^3=3, and 
5=5 is shown in Fig.I^fa). At ^ = (7r,7r), the distortion 
is a pure intracell mode, and the energy depends only 
on the intracell mode modulus B. Therefore, as shown 
in Fig.I2fb), a;(7r,7r) vanishes without the intracell mode 
(5=0), which is unphysical. 

We apply our formalism to obtain the domain wall so- 
lution for the atomic displacements between two homo- 
geneous strain states (a "twin boundary" ) due to a phase 
transition to a rectangular lattice. We then compare the 
solution to that obtained from continuum theory, where 
discreteness effects are neglected 4J . With elastic energy 

E(^l = ^-U'.e^iif + ^F^esiir, (6) 

i 

the degenerate ground state of E'rec is a uniform state 
with es=±y/ A'^/Fs^ and ei=e2=5+=5_=0. To study the 
domain wall between these two degenerate rectangular 
ground states, we consider es{i) as the order parameter 
and minimize Erll with respect to the other variables, 
using the constraint equations [Eqs. (1) and (2)] and the 
method of Lagrange multipliers. We obtain E^H = 

^^s{—k)U {k)es{k) ^ where U{k) is given in Ref. 0. 
With kx=k cos and ky=ksmO^ the expansion of 
U{k,0) about k=0 yields U{k,0) = Uo{0) + U2{0)k'^ + 
0{k^), where Uo{0) = A1A2 cos^ 2(9/(Ai sin^ 2i9 + A2), 
and U2{0) = sin^ 2i9 [6^1^25 sin^ 2(9 + 4AiA2(Ai + 
A2) cos^ 2e^3B{Al^Al sin^ 2i9)] / [24(^2 + Ai sin^ 20)^]. 
The term Uq is purely orientation-dependent without a 
length scale, and is minimized at =45^ and 135^, as ob- 
tained in Ref. . The difference between continuum and 
our discrete theories lies in the k'^ term: continuum the- 
ory commonly assumes isotropic gradients in the order 
parameter, i.e., (Vea)^ 4], whereas U2{0) is anisotropic. 
The two origins of the anisotropy are: (a) the compati- 
bility relation, Eq. (1), which has higher powers in k than 
Eq. (3) due to discreteness, and (b) the presence of shuf- 
fie mode energy. The latter can be written as gradients 
of strains, but with corrections to the phenomenological 
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FIG. 3: (Color) Atomic scale 135^ domain wall profile for 
critical length scale, Ac < 1 along the direction perpendicular 
to the domain wall: (a) strain 63 and shuffle s-, (b) differences 
in 63 (^e3=e3,atomic- 63, continuum), S- (Ss-) and displacement 
parallel to the domain wall direction {Sd\\) between the results 

from continuum theory for ^ ^ and our model that includes 
discreteness. Parameter values are Ai = 5, A2 = 4, = 4, 
B = 5, and F3 50. 

isotropic term, (Vea)^, commonly used in LG theory. As 
U2{0) is minimized for = 0^ and 90^, it competes with 
Uo{0) which prefers 6 = 45^ and 135^. Thus, the domain 
wall direction depends on the length scale with a criti- 
cal length scale Ac ~ ^/B/Al. If Xc < 1, i.e., less than 
the interatomic spacing, the domain wall has direction 
45^ or 135^ down to atomic scales. If Ac > 1, then for 
length scales smaller (larger) than Ac, the domain wall 
direction is 0^ or 90^ (45^ or 135^) and the domain wall 
has multiscale attributes. 

We examine first the case Ac < 1 that would ap- 
ply to materials with relatively large bulk modulus Ai 
('hard' materials) for fixed B. Here kx — :^ky and 
U{k) = B{1 — coskx)/{l + coskx). We illustrate the 
domain wall solution with 135^ domain wall direction. 
The only non-zero distortion modes are 63 and S- (s+ 
for a 45^ domain wall). The strain 63 reverses sign at the 
domain wall, the intracell mode S- is confined within the 
domain wall, and the atomic displacements are parallel 
to the domain wall direction. The numerical solution for 
63 and S- along a line perpendicular to the wall is shown 
in Fig. inija), for which Ac ~ 1. The e3-field and the 
corresponding displacement field near the center of the 
domain wall are shown in Fig. ^a) and^b), in which 
the red and blue colors show regions with 63 positive and 
negative, respectively. Both figures show that the cen- 
ter of the domain wall is located at bonds rather than 
sites to avoid the higher energy state of 63=0 and large 
S-. In Fig. 3(b) we compare our results with contin- 
uum theory, which predicts 63 = 63^^^ tanh(zs/^) Q and 
S- = des/2dis from Eq. (4), where is = ix^^- The dif- 
ferences in the interface region, shown in Fig. IHfb), are 
of the order of 10% of eg^^^ = ^A'^/F^,. 

The domain wall solution for Ac > 1, typical for small 
bulk modulus Ai or 'soft' materials, is shown in Fig. [^1 
for which Ac ~ a/5. The 63 field in Fig. Efa) shows that 
on length scales of the size of the system (larger than 
Ac), the diagonal orientation is still preferred. However, 




FIG. 4: Strain-e3 mode for a periodic twinned microstruc- 
ture (a) and a single defect (c) together with their correspond- 
ing displacement configurations [(b) and (d)] within the high- 
lighted window. A/"!, = 32 X 32. 




FIG. 5: (Color) Atomic scale domain wall solution for ma- 
terials with Ac > 1. Parameter values are Ai—1^ A2 — 4, 
A3 = 4, 5 = 5, and F3 = 50. Strain 62 is zero. Dark red 
corresponds to 0.28, dark blue to -0.28, and green to ~ 0. 

this diagonal domain wall consists of a 'staircase' of 0^ 
and 90^ domain walls of length scale Ac. More details 
on induced ei, s+, and S- fields around the 'staircase' 
wall in Figs.ISJb)-ISJd), their implication for functionality 
of the domain walls, and the displacement pattern are 
discussed in Ref. p]. 

A similar approach is used to find elastic texture 
around structural defects. We consider impurity atoms 
at the centers of the square of four atoms, which couples 
to the ei mode distortion of the four nearest neighbor 
atoms. The corresponding energy expression is Egq^imp = 
Esq.iat + Eimp, whcrc Esq^iat is Eq. © and Eimp is 

Eimp = ^Ciei{i)hi{i). (7) 

i 

Here, hi{i) is 1 if there is a defect at the site at 
i + (1/2,1/2), and zero otherwise. Ci represents the 
strength of the coupling. Esq^imp is minimized about ei, 
62, 63, 5+, and S- with constraints among them for given 
/ii, which gives the relations between the relaxed strain 
fields and the hi field. Explicit expressions of these re- 
lations will be presented elsewhere. As a simple case, 
we show the elastic texture around a single defect in 
Figs.Efc) andlHd). 

To illustrate the influence of lattice deformation on 
electronic properties, we couple the twin boundary and 
defect solutions obtained above with the electronic de- 
gree of freedoms in a model of superconductors. The 
electronic model Hamiltonian is defined on a square lat- 
tice: 
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(8) 



Here c^cr annihilates an electron of spin a on site i. The 
quantities and /i are the on-site impurity potential 
(if any) and the chemical potential, respectively. The 
hopping integral iij is modified by the lattice distor- 
tion. The electron-lattice coupling is approximated by 
tij = t^j [1 — acij]^ where t^j is the bare hopping integral, 
Cij is the lattice-distortion variable, and a is the coupling 
constant. In our nearest neighbor realization, the bare 
hopping integral t'^j is t for nearest neighbor sites and zero 
otherwise. Specifically, we take the form of the lattice dis- 
tortion to be: Cij = [|(R_^-+d_^-)-(Ri + di)|/|R_^--Ri|-l], 
where {Ri} are the undistorted lattice coordinates and 
{di} the lattice displacement vectors with respect to 
{Hi}. We assume an effective superconducting gap func- 
tion given by Aij = -^(QiCji — Q|Cj|), where Uij = U dij 
(i.e., attractive Hubbard-t/ model) for 5- wave supercon- 
ductivity and Uij = VSi-^^j (with 7 specifying the near- 
est neighbors to the i-th site) for (i-wave superconductiv- 
ity. By performing a Bogoliubov-Valatin transformation, 
we may diagonalize Eq. Q by solving the Bogoliubov-de 
Gennes (BdG) equation 




FIG. 6: Spatial variation of the SO OP for periodic twin 
boundaries displayed in Fig. EI a) — (a) The s-wave OP in an 
s-wave superconductor, and (b) the d-wave and (c) extended 
s-wave components of the OP in a d-wave superconductor. 
The electron- lattice coupling constant a = 3. 




(c) - 



+ 



+ 



FIG. 7: Spatial variation of the SC OP for a single defect dis- 
played in Fig. EJc) — (a) The s-wave OP in an s-wave super- 
conductor, and (b) the d-wave and (c) extended s-wave com- 
ponents of the OP in a d-wave superconductor. The electron- 
lattice coupling constant a = 3. 
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(9) 



subject to the self-consistency conditions for the super- 
conducting (SG) order parameter (OP): 



3 



<*^i^)tanh 



\2kBTj 



(10) 



Here the single particle Hamiltonian reads Hij = —iij + 
(e^ — ii)5ij. We numerically solve the BdG equations self- 
consistently. Below, we report results for two types of 
local lattice distortions at zero temperature — a super- 
lattice formed by twin boundaries and a single defect . 
We measure the length and energy in units of ao (the 
undistorted lattice constant) and t. The chemical po- 
tential fi = {) and no extrinsic impurity scattering is in- 
troduced (e^ = 0). The pairing interaction for both the 
5-wave {U) and d-wave {V) superconductors is taken to 
be 3. The typical system size is Nl = 32 x 32 with peri- 
odic boundary conditions. When the local quasiparticle 
density of states (LDOS) is computed, we implement a 
much larger system using the above small system as a 
supercell. 

In Fig. El we show the spatial variation of the SC OP 
induced by the deformation of Fig. ^a) in both s and d- 
wave superconductors. In both cases, the OP is lowered 
within the domain and is elevated at the domain wall 
(Fig. lEIa-b)). The magnitude of the OP is depressed in 
comparison to an undistorted square lattice since the lat- 
tice deformation changes the band structure, leading to a 



reduction in normal density of states at the Fermi energy. 
Even at the domain wall, where the strain induced de- 
formation is weakest, the amplitude of the enhanced OP 
is smaller than its value in an undistorted square lattice. 
This is due to the confinement from the two neighboring 
domains. In a twinned domain of a (i-wave SG, a sub- 
dominant extended s-wave component is generated in a 
real combination d ± s. Because the symmetries of two 
twinned domains are reflected into each other with re- 
spect to the twin boundary, the relative phase between 
the d- and s-wave components switches by tt when a twin 
boundary is crossed (Fig.|Hfc)). 

We show in Fig.[3the spatial variation of the supercon- 
ducting OP around the single defect (Fig. ^c)) in both 
the s-wave and d-wave superconductor cases. The OP 
is depressed at the center of the defect, and reaches its 
defect-free bulk value at the scale of the superconducting 
coherence length. Notice that for a lattice-deformation 
defect, which affects the local electron hopping integral, 
the OP has a minimum at four sites surrounding the de- 
fect center. It is different from the case of an externally 
substituted unitary impurity, where the minimum OP is 
located only at the impurity site itself The range 

of influence of such a defect can be very large depending 
on the strength of electron-lattice coupling — the elastic- 
ity propagates the electronic response. The (i-wave energy 
gap has a sign change at the nodal directions of the essen- 
tially cylindrical Fermi surface, but the d-wave OP does 
not exhibit such a sign change in real space. With the de- 
fect, an extended s-wave component of the OP is induced 
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FIG. 8: The LDOS at a twin boundary in s-wave (a) and d- 
wave (b) superconductors. Also shown are the LDOS (black 
lines) for a uniform domain. The electron-lattice coupling 
constant a = 3. 



when the dominant d-wave component is depressed at the 
defect. Strikingly, the induced 5- wave component has a 
sign change across the diagonals of the square lattice, i.e., 
sgn[cos(26>)], where 6 is the azimuthal angle with respect 
to the crystalline x axis. This is a direct manifestation 
of the d-wave pairing symmetry in real space. 

Once the self-consistency for the order parameter is 
obtained, we calculate the LDOS: 



Pi{E) 



[|<|V'(^-^J + I<|V'(^+^J] , (11) 



where f'{E) is the derivative of the Fermi distribution 
function with respect to the energy. The LDOS deter- 
mines the differential tunneling conductance, measurable 
by STM experiments 11]. Figure |H1 shows the LDOS at 
a domain wall for both types of superconductors, where 
the modulation of the OP forms a superlattice, with 
maximum at the domain wall playing the role of an off- 
diagonal potential barrier (A^j in Eq. Q). For an 5-wave 
superconductor, the quasiparticles are gapped away with 
their energy below the minimum SC OP. Outside the 
minimum of the pair potential, energy bands are formed 
by the quasiparticle scattering off the off-diagonal energy 
barriers at the domain walls. Interestingly, the bottom 
of the oscillation pattern follows the LDOS (black line) 
of a system formed by a uniform rectangular domain. 
Similar oscillations are obtained for the d-wave super- 
conductor. However, the bottom of the oscillations do 
not follow the single domain DOS (black line). In addi- 
tion, weak subgap peaks (labeled by arrows in Fig.|Hfb)) 
appear symmetrically in the LDOS on the domain wall 
but are absent in the single-domain LDOS. We speculate 
that these resonant states are due to the gradient of the 
5-wave gap component induced inside the domain. 

In Fig. El we show the calculated LDOS near the cen- 
ter of a single defect. The depression of the SC OP at 
the defect makes a quantum- well- like profile of the energy 
gap. The size and depth of the well is determined by the 
electron-lattice coupling constant. In the s-wave super- 
conductor, the well is shallow and small for weak cou- 
pling, which cannot trap low-lying quasiparticle bound 




FIG. 9: The LDOS near the center of a defect in s-wave (left 
column) and d-wave (right column) superconductors. The 
distance of the measured point away from the defect is labeled 
by its coordinate. The electron-lattice coupling constants are 
a = 3 (red lines) and 10 (blue lines). Also shown is the 
defect-free LDOS (black lines). 



states; for strong coupling constants, the well is deep 
and large so that subgap quasiparticle bound states are 
induced (the red and blue lines of Figs. Efa) and (c)). 
The energy of these low-lying states must be inbetween 
the bottom and edge of the well. Therefore, it is no- 
table that the energy of these subgap states is shifted 
toward the Fermi surface as the electron-lattice coupling 
is increased (the blue line in Figs. IHIa) and (c)). The 
electronic structure at the defect in a d-wave supercon- 
ductor becomes even richer: For a = 3 (weak coupling 
as compared to the band width of the uniform square 
lattice), the lattice distortion plays the role of a weak 
defect for the quasiparticle scattering. In this case, a res- 
onant peak with a dip exactly at the Fermi energy is seen 
(the red line in Fig.ElJb)). The overall peak comes from 
the scattering of quasiparticles off the single-particle off- 
diagonal potential (i.e., local change of the hopping in- 
tegral as a response to the lattice deformation). This 
lattice-deformation induced resonance state also exhibits 
Friedel oscillations. Typically, the peak structure ap- 
pears in the LDOS at (0,0) (We label the four sites sur- 
rounding the defect center by (0,0), (1,0), (1,1), (0,1)) 
and (-2,-2). For a = 10 (strong coupling), the 'resonant' 
peaks are pushed to higher energies (!::^ ±0.3) (the blue 
line of Fig. ini^d)). Furthermore, small shoulders appear 
close to the Fermi energy (the blue lines of Figs. Eljb) 
and (d)), which are precursors of new Andreev resonance 
states. We have also computed the LDOS without im- 
posing self-consistency on the OP and found that the 
double-peak structure is V-shaped with no existence of 
the shoulders. Therefore, the new Andreev resonance 
states originate from the confinement of the induced s- 
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wave OP. All these features are unique to an elastic defect 
in a (i-wave superconductor with short coherence length. 

In summary, we have reported an approach to "atomic- 
scale elasticity" , which uses symmetry modes of elemen- 
tary objects of atoms as distortion variables. A gradient 
expansion for the energy with anisotropic coefficients has 
been obtained, with corrections to the usual phenomeno- 
logical isotropic gradient terms used in LG theory. As an 
illustration, we obtained domain wall (twin boundary) 
solutions and elastic texture around a defect in terms of 
strain and intracell modes, and showed how the domain 
wall solutions differ from the continuum elastic soliton 
solution Using the atomic scale proffies of elastic 
texture, we studied the effects of elastic lattice deforma- 
tion on the nanoscale electronic structure in supercon- 
ductors within a BdG approach. We showed that the SC 
OP is depressed in the regions where the lattice defor- 
mation exists. The calculated LDOS suggests that the 
electronic structure is strongly modulated in response 
to the lattice deformation. In particular, it is possible 
to trap low-lying quasiparticle states around the defects. 
Images of these states will manifest the underlying long- 
range anisotropic elastic lattice deformation. These pre- 
dictions can be directly tested by STM experiments in 
new superconducting materials. Our approach is readily 
extended to other elastic textures, SC symmetries, and 
lattices, as well as coupling to other electronic models (for 
charge-transfer, charge-density-wave, spin-density-wave, 
Jahn- Teller, etc.). 
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